Structural, electronic, and magnetic properties of a ferromegnetic semiconductor: 

Co-doped Ti0 2 rutile 
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Room-temperature ferromagnetism has been recently discovered in Co-doped Ti02 rutile. Our ah 
initio density-functional theory investigations show that the substitutional Co ions incorporated into 
Ti02 rutile tend to cluster and then the neighboring interstitial sites become energetically favorable 
for Co to reside. This suggests that a Co-doped rutile containing only substitutional Co is not 
an appropriate reference bulk system in derterming the local environment of Co in polycrystalline 
(Ti,Co)02 rutile. We also find that the interstitial Co is in the low spin state and destroys the 
spin-polarization of the surrounding substitutional Co, hence reduces the average magnetic moment 
of impurity atoms. 

PACS numbers: 75.50.Pp, 61.72. Ww, 75.30Hx 
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I. INTRODUCTION 

Spintronics has emerged expeditiously as a novel tech- 
nology where both charge and spin are harnessed for new 
functionalities in devices combining standard microelec- 
tronics with spin-dependent effect. Q, 0, 0]. For their 
practical spintronic applications, a ferromagnetic semi- 
conductor is required to be not only easily integerable but 
also having a high Curie temperature. The development 
of new materials grows rapidly. Room-temperature ferro- 
magnetism has been observed, for instance, in Mn doped 
ternary compounds such as CdGeP2 H , ZnSnAs2 @ , and 
ZnGeP2 6]. A recent discovery of room-temperature fer- 
romagnetism in Co-doped Ti02 anatase by Matsumoto et 
al. using combinatorial molecular beam epitaxy (MBE) 
technique^ has motivated intensive studies on the struc- 
tural and physical properties of this material. |sl 1^. ITol ITlj 
Transparent to visible light, this material is expected to 
be of particular importance in spin-based optoelectronic 
applications such as spin- light-emission-diode and optical 
switches operating at terahertz frequency Q. 

It is interesting to note that Matsumoto et al. ob- 
served magnetic structure only in (Ti,Co)02 anatase but 
not in rutile phasej^i seemingly suggesting that the 
magnetization of the impurity atoms depends sensitively 
on the local environment surrounding them. More re- 
cently, Park et al. have successfully grown ferromagnetic 
(Ti,Co)02 rutile films by sputtering [l3|. The Curie tem- 
perature was estimated to be above 400 K for Co con- 
tent of 12%. Since the rutile phase is thermodynamically 
more stable than the anatase phase, it may have higher 
potential for technical applications. The dissimilar be- 
haviors of films grown by different methods indicate that 
the local environment of Co could have drastic effect on 
the magnetic structure of this material. Furthermore, 
by annealing the (Ti,Co)02 samples prepared by the sol- 
gel method at different temperatures, Soo et al. found 
that the local structure arround Co could remain virtu- 
ally anatase-like while the nanocrystalline materials un- 



derwent anatase-to-rutile phase transition|l4j. This phe- 
nomenon is rather unique and certainly deserves more 
detailed investigation in order to gain a microscopic un- 
derstanding of the magnetism displayed by this intriguing 
material. 

We here point out that in extracting the information of 
the local structure of Co through the X-ray absorption 
fine structure method, Soo et al. employed a refernce 
bulk system assuming all Co ions substitute for the Ti 
sites in rutile TiO^jlJ]- ^ n a recent work, we showed by 
ab initio density functional theory study that Co could 
also occupy the interstitial sites in anatase TiC>2. The 
question arise is: Is it also the case for Co-doped rutile 
Ti02? To answer this question, we have performed a 
computational investigation on the structural, electronic, 
and magnetic properties of Co-doped Ti02 rutile. We 
show that the substitutional Co ions tend to cluster, but 
this has no remarkable effect on the magnetization of 
Co. We also find that an interstitial site will become 
energetically favorable for a Co atom (in reference to bulk 
Co) when two of its neighboring substitutional sites are 
occupied by Co. The interstitial Co is in the low spin 
state (0.94/is), but it exerts strong detrimental effect 
on the magnetization of its neighboring Co (0.76fiB — * 
— 0.03/is) and therefore reduces the magnetic moment Co 
in this magnetic semiconductor. In the following, we give 
a brief outline of our theoretical procedure, a discussion 
of the results based on bonding characteristics and its 
relevance to our understanding of the magnetism in this 
Co-doped oxide. 

The calculations were done using a spin-polarized 
version of the Vienna ab initio Simulation Package 
(VASP)[T3- Fully nonlocal Vanderbilt-type ultrasoft 
pseudopotentials[16j were used to represent the electron- 
ion interaction. Exchange correlation interactions were 
described by the Perdew-Wang 1991 generalized gradi- 
ent approximation (GGA)^]. Wave functions and aug- 
mentation charges were expanded in terms of plane waves 
with an energy cutoff of 396 and 700 eV, respectively. For 
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each system, the geometry was relaxed until the atomic 
forces were smaller than 0.03 eV/A. 

The 6-atom rutile unit cell of Ti02 (see Figure la) 
is characterized by the two lattice constants a and c 
and the internal parameter u. The positions of Ti are 
(0, 0, 0) and (1/2, 1/2, 1/2); and those of O are (±u, 
±u, 0) and (1/2+it, 1/2+w, 1/2). The optimized a, 
c, and u are 4.661 A, 2.973 A, and 0.638, in good 
agreement with experiment) 181 and previous pseudopo- 
tential GGA calculations [l^. For Co-doped cases, we 
mainly dealt with a supercell containing 16 primitive 
unit cells (2a x 2a x 2c, see Figure lb) with one or two 
Ti atoms (site Bl, or, sites Bl and B2) replaced by Co 
and with/ without an interstitial Co (site A). Base on our 
knowledge of anatase (Ti,Co)02 that volume contraction 
and expansion associated with substitutional and inter- 
stitial Co have only negligible effect on the formation 
heat of (Ti,Co)0 2 [TlI, we fixed a and c at the calculated 
bulk TiC>2 values and optimized only internal freedoms 
upon Co-doping. 

To see whether substitutional Co (Co s ) cluster or not, 
we first minimized the free energy of the cell TiyCoOig 
(Co content: 12.5%). Then we doubled this cell and 
brought the two Co s together to form Tii4Co^ 1 Co^ 2 032 
(2a x 2a x 2c). We find the latter configuration is more 
stable by 0.31eV/Co, indicating that Co clustering does 
occur at this concentration. The calculated formation en- 
ergy of an interstitial Co (Co 7 ) as a neutral defect in ref- 
erence to bulk Co, £(Tii 6 Co^0 32 )-.E(Tii 6 03 2 )-£:(Co), is 
+1.75 eV, implying a sole Co 7 is highly unstable in rutile 
Ti02- This can be understood in the bond-order— band- 
strength picture. Ti:0 + Co:0 is much weaker than 
0:Ti:0, as is evident from the comparison of the forma- 
tion heat of Ti0 2 , TiO, and CoO 20]. Nonetheless, with 
the presence of Co^, the formation energy of Co 7 , i.e., 
S(Ti 15 Co| 1 Co^032)-S(Ti 15 Co| 1 3 2)-S(Co), drops to 
+0.47 eV. According to our calculations, the formation 
heat of C0O2 in a hypothetical rutile TiC>2 structure is 
—3.46 eV, slightly higher than two times of that of CoO 
(-1.59 eV). This means that Co:0 + Co:0 is only slightly 
weaker than 0:Co:0 , thus Co 7 becomes much less unsta- 
ble. Further calculations show that when both Bl and B2 
site are occupied by Co upon clustering, site A becomes 
favorable for a Co 7 with a formation energy of —0.72 eV. 
Once again, this can be explained by the higher formation 
heat of C03O4 (—9.57 eV) than two times of that of C0O2 
(—3.46 eV). Our calculations thus reveal a 2Co +C0 7 lo- 
cal structure in Co-doped Ti02 rutile. It should be noted 
that the formation of more complicated (Co s , Co 7 ) com- 
binations is also possible depending on the growth condi- 
tion, but the formation probability should be lower than 
2Co s +Co 7 . 

In order to examine the effect of clustering on mag- 
netism, we performed calculations on both intra-cell 
FM and intra-cell AFM alignments of the Co pair in a 
Tii4Co^ 1 Co|| 2 032 unit cell (For simplicity, we denote 
this cell as Tii4Co20 32 in the following) while keeping 
the inter-cell coupling ferromagnetic. The AFM phase 



is found to be unstable and it switches to FM during 
self-consistency iterations. When the two Co ions are far 
apart by 7.23 Awith one at (0, 0, 0) and the other at 
(a, a, c), the FM phase is lower in totl energy by 0.26 
eV/Co than the AFM phase. In the case of 2Co s +Co 7 , 
we checked all the three possible intra-cell spin align- 
ments of the three Co ions. All of them converge to the 
same final state, with both Co 7 and Co in the low spin 
state. The calculated total spin magnetic moment in a 
unit cell of Tii 5 Coi0 3 2, Tii4Co 2 03 2 , and Tii4Co 3 3 2, 
and that in the atomic sphere of Co (radius=1.10 A) are 
listed in Table I. It is seen that in all three cases. Co ions 
are in the low spin state, regardless to their oxidation 
states. The spin moment of Co s does not change much 
upon clustering. Co 7 , on the other hand, behaves rather 
differently. When coupled with Co s , Co 7 kills almost en- 
tirely the spin moments of Co s . As a result, the average 
spin magnetic moment of Co is reduced to 0.29 /j,b, less 
than a half of the non-Co 7 case. 

Figure 2 displays the calculated density of states for 
non-doped and various cases of Co-doped rutile. Like in 
anatase[8J, lll|, the Co states are mainly located in the 
energy-gap region and the O and Ti states are not much 
affected by Co doping. It is worth noting that with GGA, 
while Tii5Coi032 and Tii4Co2 032 just miss being insu- 
lators, the Fermi energy of their rutile counterparts just 
overpass the valence bands and fall into the energy gap. 
The presence of Co 7 does not change the insulating elec- 
tronic structure of this material although sole Co 7 ions 
would turn the material into a metallic system (DOS not 
shown), and as a consequence the magnetic interaction 
between Co ions at low temperature remains to be su- 
perexchange in the absence of carriers. 

In Figure 3, we compare the 3d DOS of Co in rutile 
with with those of Co in other Co oxides with different 
oxidation states such as Co0 2 (+4) and C03O4 (+3 and 
+2). C0O2, the end member of LiCo02, has recently 
been isolated[2l|. It seems to be of CdCl2 type with a 
monoclinic distortion, but positions of O remain to be 
determined. In this regard, we believe a study on a hy- 
pothetical C0O2 in Ti02 rutile structure is meaningful in 
illustrating the electronic structure of Co 4+ . For both 
C0O2 and C03O4, we have optimized their internal struc- 
tural parameters. We can see that although much more 
localized, d DOS of Co s in Tii5Coi032 resembles that of 
Co in C0O2 in both exchange splitting and crystal field 
splitting, suggesting they are in +4 oxidation state. In an 
octahedral field, ti g is lower than e g and thus Co 4+ will 
have a spin magnetic moment of 1/ib- As for the com- 
bination 2Co s +Co 7 , a close resemblance can be found 
between Co s (Co 7 ) in Ti0 2 rutile and Co 73 (Co A ) in 
C03O4. It thus provides a strong evidence that Co s and 
Co 7 in such a combination are in +3 and +2 oxidation 
states, respectively. The six d electrons of Co s just fill t 2g 
states and thus yield a vanishing spin magnetic moment. 
Unlike the interstial Co in anatase and Co A in C03O4, 
Co 7 in Ti0 2 rutile is in an octahedral field, therefore the 
lone electron in e g state gives a spin moment of I/is. 
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FIG. 1: Calculational unit cells of non-doped (panel a) and 
Co-doped Ti02 rutile (panel b). Large and small circles rep- 
resent O and Ti (Co), respectively. 



stitutional impurities by lowering its charge state and 
therefore plays an important role in the magnetic struc- 
ture of such material. 



FIG. 2: Total density of states for non-doped and various 
cases of Co-doped Ti02 rutile. Dotted vertical lines denote 
the Fermi energy. 

In summary, our ab initio density functional theory 
study reveals that in Co-doped rutile TiC>2, Co can oc- 
cupy interstitial sites, upon clustering of sbstitutional Co. 
As a result, a reference bulk system used to extract infor- 
mation of the local structure of Co in a polycrystallinc 
material should be one that containing both substitu- 
tional and interstitial Co. It is found that the interstitial 
Co destroys the spin-polarization of the neighboring sub- 
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FIG. 3: The calculated density of d states of Co in C0O2, 
C03O4, and various cases of Co-doped Ti02 rutile. Dotted 
vertical lines denote the Fermi energy. 
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